Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 3 March 2013 (MN style file v2.2) 

The COSMOS Density Field: A Reconstruction Using 
Both Weak Lensing and Galaxy Distributions 
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ABSTRACT 



The COSMOS field has been the subject of a wide range of observations, with a 



number of studies focusing on reconstructing the 3D dark matter density field. Typ- 
ically, these studies have focused on one given method or tracer. In this paper, we 
reconstruct the distribution of mass in the COSMOS field out to a redshift z = 1 
by combining Hubble Space Telescope weak lensing measurements with zCOSMOS 
spectroscopic measurements of galaxy clustering. The distribution of galaxies traces 
the distribution of mass with high resolution (particularly in redshift, which is not 
possible with lensing), and the lensing data empirically calibrates the mass normali- 
sation (bypassing the need for theoretical models) . Two steps are needed to convert a 
galaxy survey into a density field. The first step is to create a smooth field from the 
galaxy positions, which is a point field. We investigate four possible methods for this: 
(i) Gaussian smoothing, (ii) convolution with truncated isothermal sphere, (iii) fifth 
nearest neighbour smoothing and (iv) a muliti-scale entropy method. The second step 
is to rescale this density field using a bias prescription. We calculate the optimal bias 
scaling for each method by comparing predictions from the smoothed density field 
with the measured weak lensing data, on a galaxy-by-galaxy basis. In general, we find 
scale-independent bias for all the smoothing schemes, to a precision of 10%. For the 
nearest neighbour smoothing case, we find the bias to be 2.51 ± 0.25. We also find 
evidence for a strongly evolving bias, increasing by a factor of ~ 3.5 between redshifts 
< z < 0.8. We believe this strong evolution can be explained by the fact that we use 
a flux limited sample to build the density field. 
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1 INTRODUCTION 



The COSMOS field (|Scoville et all 120071 ) is the largest 
region of the sky that has been mapped contiguously 
with the Hubble Space Telescope (HST). This field has 
since been the subject of a wide range of studies aimed 
at measuring the detailed properties of objects in this 
field in m any ways, including spe ctroscopic follow-up with 
the VL T dLillv et all |2007| . l2009h . infrared imaging with 
Spitzer (ISanders et al.ll2007), x-ray observations with X MM 



and Chandara (Elvis et al 



20071; lHasinger et al.ll2007n , UV 



imaging with GALEX (IZamoiski et al.ll2007l )T Subaru imag- 



ing (ITaniguchi et al]|2007l) an d radio observations with the 
VLA (jSchinnerer et al.l 12007). The wealth of data coming 
from this region makes the COSMOS field perfect for de- 
veloping techniques that bring together different data sets. 
These probe combination methods are useful for current 
studies, since it allows us to maximise the information com- 
ing from current data, as well as helping us to prepare for 
future wide-fi eld surveys that will increasingly rely on probe 
combination i|Albrecht et all 120061 ; IPeacock et~aLl 120061 ) to 
make high precision measurements. For instance such sur- 
veys will aim to measure the dark energy properties, such 
as the equation of state, at the percent level. 

We will focus in this paper on methods for reconstruct- 
ing the density field from weak lensing and the spatial distri- 
bution of galaxies. This type of density field reconstruction 
work is a very active field an d has been performed on a num- 
ber of surveys. For instance. iKitaura et al.l (|2009l ) performed 
a density reconstruction of the Sloan Digital Sky Survey 
(SPS S) north cap using data release 6 (|Adelman- McCarthy! 
2008) using Wiener Filtering. A number of density field re- 
construction studies have also been performed on the COS- 
MOS field. These include reco nstructions of the density field 
using the gal axy distribution (IKovac et al. l l2010l . l201ll ) and 
weak lensing (M assev et al. I l2007l ). Each of these approaches 
has its own strengths and drawbacks. For instance, the den- 
sity fields constructed from the weak lensing data alone have 
very poor resolution in redshift. This is due to a convolu- 
tion along the line of sight by 'the lensing efficiency func- 
tion', which we discuss in more detail in section 12.21 The 
strength of the lensing maps, however, is that the weak lens- 
ing signal is a direct probe of the underlying matter. Maps 
constructed using galaxy positions rely on the galaxies act- 
ing as tracers of the matter field (see section 12. f \ . On the 
positive side, these galaxy position reconstructions have sub- 
stantially better resolution in redshift than what is possible 
with weak lensing. Nonetheless, since the galaxies that we 
see are, at best, biased tracers of the underlying density 
field, additional assumptions and simplifications are needed 
to produce a reconstruction of the density field. 

The work presented here, therefore, aims to perform a 
matter reconstruction using both the galaxy position and 
weak lensing data. Such a combination allows us to con- 
struct a density field with high resolution in redshift that 
is calibrated from the data. We do this by measuring the 
expected lensing signal for each galaxy in the weak lensing 
survey that would come from a particular density field re- 
construction. With this we can then look for correlations 
between the measured shear signals and predictions on a 
galaxy by galaxy basis. This approach to the problem is 
very powerful because it gives us freedom to choose how we 



average the data to reduce noise. For instance, we do not 
need to bin the data spatially to average out intrinsic shape 
noise. Instead we can average over predicted shear, which is 
both easier and more stable, or we can do away with binning 
entirely through direct correlations. 

This paper is organised as follows. In section [2] we 
present an overview of some of the main issues associated 
with density reconstruction for galaxy clustering and weak 
lensing. We then present the COSMOS data sets in section 
[3] with a discussion of our methodology in section 3] We give 
our results and conclusions in sections [5] and [6] 



2 OVERVIEW OF THE ISSUES 

When studying the density field, p, we typically focus our 
attention, in cosmology, on the over-density, 8. For a given 
cosmic time t this is dened as 



8 = 



— 1 

P 



(1) 



where p is the mean cosmic density at that time. We do this 
because it is these perturbations of the smooth background 
density that drive structure formation and cause the bending 
of light rays. 



2.1 Galaxy Bias 

It is well known that since galaxies form at peaks in the back- 
ground density field, the statistics of their distribution follow 
that of the underlying dark matter in a biased way. For in- 
stance, the two-point correlation of the peaks is boosted rel- 
ative to the two-point correlation function of the underlying 
field, where the boost factor depends on the threshold used 
to define a peak. This boosting factor in the correlation func- 
tio n, which we w ill call Kaiser bias (6fc), was first identified 
bv lKaiserll|l984 ). but it has since been calcu l ated by a num- 
ber of other works jEfstathiou et alj 
1 19891 : lKaiserill987l : IPeacock fc Heavens! 



1988 



ICole fc Kaise: 



1985; Bardccn 198' 



for Gaussian random fields, as well as bein g studied us- 
ing pr ocess named peak-background splitting ( Manera et al.l 
l|2010l). Since, for instanc e, in the Press-Schechter framework 
l|Press fc Schechterill974l ). the peaks above a critical thresh- 
old become halos and galaxies, the bias is often used to refer 
to the factor that links the galaxy and dark matter correla- 
tion functions, 



Mr) = b 2 k (r)U(r), 



(2) 



where £ 9 and £ m are the two-point correlation functions of 
the galaxies and matter, respectively, as a function of sepa- 
ration distance r at a given redshift. In principle, the Kaiser 
bias, fefe, could be scale dependent. With the bias defined 
in this way, its meaning and interpretation are straightfor- 
ward and unambiguous. Both correlations can be measured 
in numerical simulations (assuming that the halos host the 
galaxies) . This has been done by many autho rs (for exam- 
ples, see lWeinberg et al.|[20oi ; iMo et al.lll996l) . The general 
findings are that on large scales the bias is constant and re- 
sults are consistent with expectation s from the Halo model 
l|Mo fc White! 1 19961 : iMa fc Frvl I2OO0I ) . which links the halo 
mass to density threshold of the matter fields. 
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The link between the galaxy population and the under- 
lying density field in real space is more complex due to the 
fact that galaxies form a set of points and not a smooth field. 
We then need to translate a set of coordinates in space (for 
each galaxy we have redshift z and the two coordinates 6 
and <f>) into a real space galaxy density measure (S 3 ), which 
can then be used to build the smooth underlying matter 
density field (<5 m ). The steps, therefore are, 



{9 t , fa, zi) -> S 9 (6, <t>, z) -> S m (6, <t>, z). 



(3) 



The relationship between the smoothed galaxy density field 
and the underlying matter field will then depend on the 
way that the smooth galaxy field was created. This means 
that for real space density reconstruction there will not be 
a unique bias. Instead the bias (&x), where we use X to de- 
note the smoothing method, will depend on the particular 
scheme that has been chosen for going from galaxy positions 
to a continuous density distribution. For instance, produc- 
ing a continuous galaxy field by smoothing with a Gaussian 
of fixed angular size, i.e. convolving a Gaussian and a set of 
delta functions, will likely require a different scaling (&Gauss) 
than a scheme with a dynamic smoothing scale, such as es- 
timating the density using a fixed number of nearest neigh- 
bours (&nn). 

The simplest link between the continuous galaxy den- 
sity field and the matter over-density is linear bias, where 

8 9 x =b x S m . (4) 

This simple relation is similar to the Kaiser bias in equa- 
tion |2]l. except that the small scales will be affected by 
the smoothing used to go from galaxy points to the smooth 
density field. It is also possible to invoke more complex rela- 
tionships, such as a non-linear and redshift dependent bias, 



6i(2)<5 m + b 2 (z)(S m f 



(•») 



and in general the bias can be expected to be redshift depen- 
dent. Note that although we have omitted the subscript X, 
which would cause our notation to become clumsy, the bi- 
ases are still specific to a given smoothing scheme. This will 
be true for all real space biases (and inverse biases) used in 
this paper even if we omit the X subscript for convenience. 

The nonlinear bias in equation [5] would likely lead to a 
scale dependent Kaiser bias. In all cases, the most accurate 
bias in real space will depend on (i) the tracers used, (ii) 
the smoothing scheme for going between galaxies as points 
and a smooth density field and (iii) the way that the density 
field will be used. On the last point, it is not clear that there 
exists a single density field that is optimal for all users. A 
density field to be used for cosmology may have different re- 
quirements than one used for galaxy evolution studies. For 
instance, for the former it may be best to use a simple fil- 
ter so that the statistics of the matter field can easily be 
compared with prediction from theory, whereas the complex 
mapping from point to smooth field, e.g. involving dynamic 
smoothing scales, may be more optimal for the latter. 

2.2 Weak Lensing 

Gravitational lensing effects are caused by the bending of 
light by intervening matter as it travels from source to ob- 
server. Because of this, lensing observables are integrated 
quantities along the line of sight. It is then difficult to make 



accurate measurements in the redshift direction. The ad- 
vantage of weak lensing, however, is that the signal does not 
depend on the type of matter along the light's trajectory. 
This means that dark matter, which is the dominant matter 
component, dominates the lensing signal. For this reason, 
lensing reconstructions of the matter den sity have focused 
mainly on 2D projected reconstructions I Kaiser fe Sq uires 



Soucail20q7l;lBerge et al.ll20oijlkubo et al 

2008: 



1993; Gavazzi 

20091 : iHevmans et all l200Sl : iKurtz et alJ 120111 ). Although 
the redshift dependence is weak, there have also been at- 
tempts to perform fu ll 3D reconstructions usi n g lensing data 
llHu fc Keetonl 120(321: bacon fc Tavlorl 120031 : Ikassev et all 
120071 ; ISimon et alj|200a ). These methods are based on weak 
lensing tomography, which employs the sensitivity of weak 
lensing signal to the source redshift. 

In weak lensing, the most commonly used observable 
is cosmic shear. This is where the bending of the light rays 
causes the images of background galaxies to become sheared. 
For instance, if a background galaxy is initially a circle, cos- 
mic shear would distort the image into an ellipse. The effect 
is very subtle (causing typically percent level changes in the 
axis ratios of galaxies), so it must be studied statistically 
over very large samples of galaxies. Here, we give a very brief 
summary of the lensing basics and point the reader to th e 
following revi e w art i cles for more details: iRefregierl &2 003): 
iMunshi et ail l^OOSMBartelmann fc Schneider! ^OOll ). 

The lensing induced effect that warps a galaxy image 
can be described to first order using the distortion matrix A 



A = (1-k) 



-91 

'92 



-92 



(6) 



where k, is known as the convergence, since the first term 
causes a change in the size of the image, and g\ and <?2 are 
the two components of what is called the reduced cosmic 
shear, which causes an anisotropic shearing of the galaxy 
image. The reduced shear is related to the shear 7 through 
the expression 



9 = 



7 



1 - ft 



(7) 



Often in weak lensing we make the approximation that 
7 = g. However, it is worth noting that the actual observ- 
able is always reduced shear. The method we present here 
is able to properly account for the reduced shear. We find 
that for the accuracy possible with the COSMOS survey, 
approximating the reduced shear with the shear works well, 
but the distinction is likely to become more important for 
future surveys. 

Both shear and convergence can be related to each other 
through their dependence on the lensing potential ip, 

« = (V'.n + V>,22)/ 2 ; 71 = C0,n - '0,22)/2; 72 = ip,i2, (8) 

where the subscripts denote a second order partial derivative 
such that ij) t ij — d 2 ip/99id8j. Finally, in the weak limit typ- 
ical of cosmic shear studies, where the light path is mildly 
perturbed, the convergence can be linked directly to the 
mass through an integration along the line of sight, 



2c 2 



x(Xo -X) s 



dx, 



(9) 



where x i s the comoving radial distance and Xs i s the dis- 
tance to the source. This relies on a simplification known as 
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the Born approximation. We can see from this expression 
that the strength of the lensing signal from a given over- 
density depends on the redshift of the background source 
that is being lensed. In this way, the lensing signal can be 
seen as a radial convolution of the density with the lensing 
efficiency function. Methods have been put forward for re- 
covering this radial information and eff ectively performing 
a deconvolution fsee lSimon et al.l (2009) for further discus- 
sion). However, as with any deconvolution of noisy data, a 
perfect recovery of the original density field is not possible 
and the weak lensing reconstructed maps have very low res- 
olution in redshift. 



2.3 Joint Analysis 

There are a number of ways that information coming from 
weak lensing data and the spatial distributions of galaxies 
can be combined. Two examples that are worth outlining 
are ones based on comparing statistical properties, such 
as correlation functions computed for galaxies and for the 
shear, and what is known as galaxy-galaxy lensing, which 
depends on the cross-correlation between the position of a 
foreground galaxy and the lensing of a bac k ground galaxy 
dBrainerd et all (l99rj; iMcKav et alj l200ll; Isheldon et al.l 



2004 iMandelbaum et all 120061; IJohnston et all 120071 ; 
Sheldon et al l 120091 ; ILeauthaud et al l I2010T) . The former 
approach is most useful when studying the galaxy biasing 
using a Kaiser-like bias (equation (2j. This results natu- 
rally by measuring the correlation function of the galaxy 
distribution and comparing it to the correlation of the 
lensing, both of which can be measured directly. It is also 
useful in such an approach to compare different types of 
correlation functions, such as the aperture mass statistic 



llHoekstra et al.ll2002l ; |Schneidedli998r i. This is the approach 
taken by I Julio et al.1 (f2012T ) in studying the bias properties 
of galaxies in the COSMOS field. 

The second approach of using galaxy-galaxy lensing has 
also been widely studied. An exam ple of galaxy-gala x y lens - 
ing in COSMOS is presented by ILeauthaud et al.1 (201ll ). 
The strength of such a study is its ability to measure the 
average properties, such as the radial profile, of a given set 
of galaxies. However, the lensing effect of known secondary 
lenses along the line of sight typically is not included in 
the data analysis, instead it is included statistically as a 
'2-HALO term'. 

Our aim here is to create a real space density recon- 
struction. We take a more direct approach, outlined below 
in section |3J to relate and calibrate the density of galax- 
ies and the overall matter density field in real space. Our 
method, therefore, sits between the pure correlation func- 
tion approach and the classic galaxy-galaxy lensing studies. 



3 COMOS/zCOMOS DATA 
3.1 The COSMOS Field 

Much of the data from the COSMOS field that we use here 
has already been presented and described in detail in the 
literature. Therefore, in the sections that follow we highlight 
some of the key features of this data, and we refer the reader 
to the appropriate sources for more of the technical details. 



3.2 Weak Lensing Data: Hubble ACS Imaging 

Our weak lensing analysis relies on detailed measurements 
of the shapes of the galaxies in the COSMOS field. For this 
we use the shape catalogs generated from t he HST ACS im- 
ages. T he ca talogs are desc ri bed in detail in lLeauthaud et all 
(2007) and iRhodes et~aH (|2007T ). with the pipeline hav- 
ing received a number of improvements since the n. These 
updates are described in ILeauthaud et all (|201 lh and in- 
clude improvement s in the treatment of Charge Transfer 
Inefficiency (CTI) (jMassev et alj|20ich . The lensing cata- 
log is constructed from 575 ACS/WFC tiles with a total of 
1.2 x 10 6 galaxies to a limiting magnitude of If&uw = 26.5. 
After lensing selection cuts are made, the final COSMOS 
weak lensing catalog contains 3.9 x 10 5 galaxies. These have 
accurate shape measurements and correspond to a density 
of 66 galaxies per square arc minute over the 1.64 square 
degree area covered by the lensing catalog. Figure [1] shows 
the redshift distribution of our lensed galaxies along with 
curves showing the lensing efficiency function, see equations 
[9] and [10] for each of the samples. 



3.3 Weak Lensing Data: Photometric Redshifts 

The redshifts of the lensed galaxy sample have been mea- 
sured through their photometry. For the work p resented 
here, we use v 1.8 (updated from lllbert et al.1 120091 '). These 
redshifts were determined using the 30 band mult iwave- 
length data analysis presented in lllbert et al. (2009). This 
included deep Ks, J, and u band data, which allows for ac- 
curate photo-z measurements at z > 1 through the 4000A 
break. Mor e details on t he da ta and the photometry can 
be found in ICapak et al] l|2007h . The photo-z measurements 
used a template fitting method (Le Phare) that was cali- 
brated with large spe ctroscopic sample s from VLT-VIMO S 
and Keck-DEIMOS l|Lillv et al.l 120071 ; lllbert et al.l I2009T I. 
The dispersion in the photo-zs as measured by comparing 
to the spectroscopic redshifts is Az/(1 + z spec ) = 0.007 at 
iAB < 22.5, where Az = z spec — z p hot- 



3.4 Spectroscopic Redshifts 

For building the density field from the gal axy distribution , 
we re ly on the zCOSMOS galaxy sample (|Lillv et al.ll2007l . 
l2009h . We use galaxies with IAB < 22.5 in the 1.1 deg2 
region of the COSMOS field. Roughly 60% of these galax- 
ies have spectroscopically confirmed redshifts from the 20k 
zCOSMOS-bright sample, with the redshifts of the remain- 
ing galaxies having been determined photometrically and 
have typical errors of roughly 8z — 0.01(1 + z). This pro- 
vides us with a flux limited population. We have decided to 
use a flux limited sample because this has the advantage of 
containing the largest possible number of galaxies. An al- 
ternative option is to use a volume limited sample, which 
would likely have a simpler bias prescription. If the bias is 
introduced through a theory prior, this would be very use- 
ful. However, since our primary objective is to reconstruct 
the density field and we are able to measure the bias em- 
pirically using lensing, we are free to use any galaxy sample 
even though the bias may be complex. Since our measure- 
ment of the bias is a means to an end, it does not need to 
be easy to interpret. 
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Photo-z 



Figure 1. The redshift distribution of the background lensed 
galaxies. The red sample shows the galaxies with photometric 
redshift of less than one and the blue sample are the galaxies 
with photometric redshifts between one and two. The dotted, 
dashed and solid curves show the lensing efficiency function for 
the red, blue and both samples, respectively. Each of these curves 
has been normalised so that the area is one. The means of each 
curve give an effective redshift that is probed by that sample. 
These are z e // = 0.36 (red sample), z e // = 0.70 (blue sample) 



andz. 



0.62 (all galaxies). 



4 OUR METHODOLOGY 

4.1 Building a 3D Density and Lensing Cube 

Our objective in this study is to build the matter density 
field ou t to a redshif t of z = 1. This is similar to what was 
done in iKovac et af] (|20ld ). We focus on the central region 
of the COSMOS field (149.575 < R.A. < 150.675 deg and 
1.75 < DEC < 2.7 deg), which is covered by both the weak 
lensing ACS catalog and the zCOSMOS data. We divide this 
volume into a 3D grid with a resolution of 256 x 256 x 500. We 
then place galaxies into this volume based on their angular 
position and redshift. In the redshift direction we account 
for uncertainties in the redshift estimation by using the full 
probability distribution function (PDF) coming from the 
redshift estimation code. For the zCOSMOS galaxies, the 
redshift PDF is effectively a delta function, so these galax- 
ies are assigned to a single pixel, while the photo-z galaxies 
are distributed over several pixels in redshift (depending on 
the photo-z error of a given galaxy). With this 3D galaxy 
density field, we can move towards calculating the conver- 
gence and shear inside this volume. To do this, we first need 
to convert the 3D galaxy density fiel d (p 9 ) into a galax y 
over-density field (S 9 - see equation [1} (|Kovac et al.ll201lT ). 

We recall, from equation^ that the convergence is given 
by the integral of the matter over density (5™ 1 ), 



8 m (z)w(z)dz, 



(10) 



tains all the weight functions from equation [9] This can be 
converted into a predicted 3D shear 7 P field through the re- 
lations hips shown i n equ ation [5] This is similar to methods 
used in lPires et aj ([2009 ). To link this calculation with our 
galaxy over density field (S 9 ), we need only to introduce the 
relationship between this and the matter over density (8 m ), 
which we do through the bias (as discussed in section [2~TT) . 
Since the lensing observable comes from the matter over- 
density, it is simpler for our applications to reverse the usual 
bias relationship (equation [4|. Instead we will use fi = 1/6 
and work with 



5 m = fiS 9 

for constant bias and 
*"(*) = fi(z)S 9 (z) 



(11) 



(12) 



for the analog of redshift dependent bias. Once again, note 
that we have dropped the subscript X to simplify our nota- 
tion. It is also possible to work with an analog of non-linear 
bias (equation O, S m = fii(z)S 9 + fi2{z)(5 3 ) 2 . However, for 
the work presented here we focus exclusively on linear bias. 
Clearly, in the simplest case of linear biasing, then b = l//i. 
Here, we focus on reconstructions to first order (i.e. linear 
and only including S 9 terms). From this, we see that pre- 
dicted convergence is thus given by 



(13) 



re p = / fi(z)8 9 (z)w(z)dz, 

and for the case where fi is independent of redshift 



k p — fj, S 9 (z)w(z)dz = \1K 9 



(14) 



where re 9 is the convergence of the raw galaxy density field 
without a bias. We find that a parameterised expansion of 
H in terms of z/(l + z), such that /i(z) = fio + [s.\zj(\ + z), 
gives a convenient form to explore redshift evolution. In this 
case, the convergence can be separated as 



P 9 i f 9 

K = /ioft + HlK , 

where re' 9 is defined as 



S 9 (z)w(z)dz. 



1 + 2 



(15) 



(16) 



where re p would be our predicted convergence and w(z) con- 



In this separable way, we can easily investigate the different 
contributions from a constant bias term and test whether 
an additional evolving term improves the agreement with 
the lensing data. Also, it is useful to remember at this point 
that a given convergence field can be converted into a cor- 
responding shear field using the relation shown in equation 
O In this case, k v leads to y p , re 9 — > 7 9 and re' 9 -> ~f' 9 . 



4.2 Background Cosmology 

To calculate the expected lensing signal from a given density 
field, we need to assume a background cosmology. Specifi- 
cally, we need to assume a background expansion that re- 
lates redshift and radial distances. The rationale that we 
have chosen to adopt is to separate cosmology constraints 
coming from geometry and structure growth and to only 
use data coming from geometry measures to set our fidu- 
cial cosmology. We do this because constructing a den- 
sity field primarily focuses on growth, so we try to limit 
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the cross-talk between external data and internal measure- 
ments coming from COSMOS. We have chosen our fidu- 
cial cosmology based on the analysis by the WMAP team 
(available on the NASA WMAP websit(0), which com- 
bines (i) the seven ye ar WMAP Cosmic M icrowave Back- 
ground (CMB) data ijKomatsu et al.ll201lf ). (ii) the Type 
la supernovae comp ilation from the extended SDSS dataset 
jKessler et alj|2009[) and fiii ) Baryonic Acoustic Oscillations 
(BAO) (jPercival et all201Cl ). Our fiducial cosmology param- 
eters are O m = 0.278, fi A = 0.722 and h = 0.699. The rest 
of the cosmology parameters, such as the equation of state 
to, are taken to be consistent with standard ACDM (e.g. 
io = -1). 



4.3 Investigation of the Impact of Smoothing 



The procedure that we outlined in section 14711 for converting 
the galaxies positions into a 3D grid can be seen as a form of 
smoothing. However, this is a very minimal level of smooth- 
ing. We can see this when we consider that the number of 
pixels in our grid is 256 x 256 x 500 fa 3.3 x 10 7 , which is 
close to three orders of magnitude larger than the number 
of galaxies that we place in the grid. At this level the grid 
is extremely sparse (mostly zero) and it is very difficult to 
construct convergence and shear fields (for some discussions 
of im portance of smoothing in calculating lensing properties 



ng 1 

seelAmara et alfeOOfj ; lAubert et"alll2007i ; iMetcalf fc Amaral 
120121 ). Further smoothing is therefore needed. 

We have explored four procedures for converting the 
galaxy point field into a smoothed continuous field. Each of 
these methods is a perfectly valid method for constructing 
a smooth density field from a point distribution of galax- 
ies. In this paper we do not address the question of which 
method is best since this is likely to strongly depend on the 
reason for wanting to create a density field. Instead we fo- 
cus on how lensing can be used to calibrate the density field 
once a given smoothing scheme has been chosen. Since each 
different smoothing scheme will produce a different realisa- 
tion of the smooth density field Sx , where X would denote 
the smoothing scheme, it is possible that different smooth- 
ing methods would have different optimal bias scaling (bx)- 
These smoothing methods we use are 

• Gaussian Smoothing: In this approach, the 2D mass 
field is convolved with a 2D Gaussian smoothing kernel. The 
results we show use a Gaussian with a standard deviation of 
1.2 arcmin, which corresponds to 5 pixels in the tangential 
direction. 

• Truncated Isothermal Sphere: A natural exten- 
sion to the generic Gaussian smoothing is to modify the 
smoothing kernel using the profile of a Truncated Singular 
Isothermal Sphere (TSIS). This kernel is given by £(#) = 
a 2 / 9 tan" 1 Of} - 6 2 /6), where E is the 2D mass, a is the 
velocity dispersion for a Singular Isothermal Sphere (SIS) 
and is set by the mass of the TSIS. Again, for this study, 
all TSIS's have been normalised to have an area of one. Fi- 
nally, we have set the truncation radius, 0t, to 4.7 arcmin 
(20 pixels). 

1 http: / /lambda. gsfc.nasa.gov / product /map /current / 
params /lcdm_sz_lens_wmap7_bao_snsalt . cfm 



• Neares t-Neighbor: This is the approach adopted by 
iKovac et af] (|2010h . This method works by calculating the 
density of a given point based on the distance to the fifth 
nearest neighbors. In this way, the smoothing scale is adap- 
tive and smooths over large scales in low de nsity regions 
and sm all length scales in high density regions. IKovac et al.l 
|2010l ) have also implemented schemes for edge correction 
and treatment of mask, which they discuss in detail in their 
paper. Unlike the other three smoothing schemes that we 
have explored, this NN does not use the gridded galaxy dis- 
tribution discussed in section | 4.1l I n stead we use exactly the 
same procedure as in IKovac et all (|201Gh . which we remap 
to a 256 x 256 x 500 grid so as to match the over-densities 
from our other methods. 

• Multiscale Entropy Filter: For a further alternative 
filtering scheme we also investigate the impact of using a 
multiscale entropy filter (MEF). This method has been im- 
plem ented in the software package MRLens (|Starck et al.l 
2006), which is publicly available. For those familiar with 
this method, we note that we decompose a given 2D mass 
sheet into eight wavelet scales and remove the first two scales 
(the first scale corresponds to the pixel scale). We note here 
that we have not attempted to modify the MEF implemen- 
tation to optimize it for the noise properties of our kappa 
maps. Instead, we use the routines as they are from the MR- 
Lens package. This is likely to mean that our MEF smooth- 
ing is suboptimal. We, therefore, use the MEF filter as a 
point of comparison, but we will not be able to make a gen- 
eral statement about the absolute merit of MEF-like meth- 
ods compared to our other smoothing schemes. 

In principle, each of these filters could be applied to 
either the 2D slices of the density field or the convergence 
field at each redshift, the latter being an integral quantity of 
the former. Since the aim of this present work is to construct 
the density field, we have focused on smoothing the 2D mass 
sheets. 



4.4 Using Predicted Shears 

By implementing the above approach we are able to calcu- 
late the predicted cosmic shear from a given density field 
at the position of every galaxy in the lensing catalog. For 
each of these galaxies, we would then have measured shear 
71 and 72, noting that shear has two components and that 
the expected shear coming from the density field for each 
galaxy is 7^ and 7^. Let us assume then that 

H = % + 7» , 



(17) 



where 7^ is a term that contains the noise contributions to 
the measured shear and i can take the value of 1 or 2 for 
the two shear components. The predicted shear is linked to 
the shear coming from the galaxy density field (7?) through 
the inverse bias, 7? = fj.'yf. Finally, we can extend the cal- 
culation of the expected lensing signal out to galaxies at 
redshifts higher than z = 1, i.e. outside our density data 
cube. These galaxies will be subjected to an additional lens- 
ing from structures between z = 1 and the galaxy. However, 
if this structure is not correlated to the structure inside our 
cube, then this extra lensing signal is effectively an addi- 
tional source of noise. This will likely be subdominant to 
other sources of noise, hence we ignore it. 
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4-4-1 Fits to Predicted Shear 

The average measured shear should go to zero for large num- 
ber of galaxies (i.e. {7;) — > 0). However, if we bin the mea- 
sured shear according to the predicted shear from the den- 
sity field, then ( / y)- yp —¥ j p . Therefore, for a good density 
reconstruction, a plot of (7) 7p vs. j p should, for the right 
bias, give a straight line with a gradient of one. In terms of 
X 2 , we can use the quality of the fit as 



2 

X = 



E 



((7), 



■7pJ 



<7 2 > 



(18) 



■yp 



where (7 2 ) 7p is the variance of the data for fixed 7 P bin, 
which will be dominated by the errors. This simple view is 
clear to understand if indeed the true density field is well 
described by the constant linear bias model in equation [4] 
In the case of more complex relations it is important to 
make the distinction between variations in the underlying 
density that are not correlated with the density field, which 
will be washed away in the averaging process, and deviations 
that are correlated with the galaxy density. An example of 
the latter would be a non-linear bias, such as the second 
term on the right-hand side of equation [5] For this reason, 
a more precise way to understand our procedure is that we 
are measuring the best-fit bias, using a linear approximation, 
that links the galaxy over-density field to the matter over- 
density field. 

4-4-8 Zero-lag Covariance for Constant Bias 

As well as the simple fitting procedure outlined in the pre- 
ceding section, it is also possible to make the bias measure- 
ments without needing to bin the data. We will show in the 
results section that both methods give consistent results. 
However, the advantage of removing the binning step is that 
calculations become more stable. This becomes especially 
important in the case of evolving bias, which we will discuss 
in the next section. The 'no binning' approach is to look at 
the covariances between the measured and predicted data. 
Since we have two components and two measures of shear, 
each galaxy gives a four element data vector {7® , 7|, 71, 72}. 
We manipulate the elements of the data vector such that 

(19) 



7i 



7i 



7? 



where y/ ((7? ) 2 ) is the standard deviation of each of the 
shear components over all galaxies, to create a new data 
vector for each galaxy {7* , -yf, 71, 72}. For constant fj,, the 
(zero-lag) covariance matrix of these elements over all galax- 
ies is 



/ <7?7?> <7 2 9 7?> 
- <72 9 7 2 9 > 



/' 



(71 7f) 
(717!) 
(7171 > 



/' 


1 + <72 V 72 V > 



(727?} 

<7272 9 > 

(7271) 
(7272) 



+ N, 



(20) 



(21) 



where the matrix shows the leading order terms and N is a 
noise matrix that tends to zero as the total number of galax- 



ies is increased. The covariance matrix, therefore, contains 
a wealth of information that we can use to measure \i and 
estimate its errors. 



4-4-3 Zero-lag Covariance for Evolving Bias 

For the case with an evolving galaxy bias, we can analyze 
the data in similar ways. When working with the data, we 
found that while the \ 2 method gives stable results for the 
case of constant bias, the results become unstable when the 
redshift bias is allowed to evolve. We find that this comes 
from the binning step and is due to the fact that the under- 
lying signal is much weaker than the noise per galaxy. For 
the constant bias case, the relative ranking of the predicted 
shears of different lensed galaxies is maintained as we vary 
the bias. Therefore, when we bin in predicted shear, a given 
bin will contain the same lensed galaxies. This makes the 
calculations stable when varying the bias, but it does mean 
that we need to be cautious about the overall normalization 
of the x 2 functions. In the case of the evolving bias, the rank 
order of the predicted shears changes, which can change the 
compositions of the bins and cause erratic results becausee 
of the large noise terms. Instead, we have found that an 
extension of the zero-lag covariance method gives stable re- 
sults. We also performed simple Monte-Carlo realsations of 
the data to confirm this effect. 

We can extend the zero-lag covariance method to evolv- 
ing bias case by noting that, from equation ll5l the predicted 
shear can be decomposed into two parts, 



fJ.o% +Afi7i • 



(22) 



Here 7? are the two components of the predicted shears, jj-o 
and fi\ are the parameter expansions of the inverse bias. The 
shear prediction 7® and 7^ 9 can be calculated from the den- 
sity field directly and thus can be compared directly to the 
measures shear 7^. By measuring the covariances of these 
three quantities we can construct /xo and pi, which are re- 
lated by 



Mo 



and 



(7i7 9 ) (7- 9 7f) 
<7 ! 9 7 ! 9 ) 1 <7 ! 9 7 ! 9 ) 



<7 ! ' 9 7 I ' 9 > /i °(7i 9 7i 9 )' 
We see from equations[23]and[24]that we have two equations 
with two variables, so we can easily solve for fio and (ii. The 
other advantage of this approach is that the errors on /io and 
Hi can directly be calculated from combinations of fourth 
and second order moments of the data. For a d iscussion of 
these techniques see Chapter 5 of iLuptonl (|l993l ). 



<7^ 9 7 9 > 



(23) 



(24) 



5 RESULTS 

5.1 Constant Bias 

5.1.1 Constant Bias with z < 1 

In our first test, we create a self consistent region out to 
redshift of z = 1. Inside this 3D volume we create the num- 
ber density field, using the methods outlined in section 11731 
and the shear field by integrating the gradient of the lensing 
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Figure 2. The predicted shears vs measured shears for the 
Nearest-Neighbour (NN) smoothing cases. The top panels show 
results using lensed galaxies out to z=l with a bias of 1.59, and 
the bottom panels show results that include lensed galaxies out 
to z=2, where we have set the bias to 2.51. The results for the 
two shear components are shown separately (71 on the left and 72 
on the right). To build these results we bin the data in predicted 
shear. In the upper panel, each data point is an average of ~ 8000 
galaxies, and in the lower panel each data point is averaged over 
~ 15000 galaxies. The red results show what happens when a bias 
of one is used to make the shear predictions. 



potential along the line of sight. By calculating the shear 
field inside this full volume, we are able to assign a shear to 
each of the galaxies in the COSMOS lensing catalog. The 
top panels of [2] show a comparison between the predicted 
shear, which include the best fit bias, coming from the NN 
smoothed density field and the measured shear for lensed 
galaxies with redshifts less than one. In this figure, we have 
binned the data in predicted shear so that each data point 
contains roughly eight thousand galaxies. In this way, the 
random scatter coming from intrinsic shape noise is reduced. 
We see a clear correlation between the predicted shear and 
the measured shear once the best fit bias of 1.59 is used. 
The red symbols in background show the results when a 
bias of one is used. Note that the two shear components 
give separate measures of the lensing signal. It is therefore 
very important that both give consistent results in figure 
[2] This is equivalent to noting that the covariance matrix 
in equation [21] contains two independent entries for /1, one 
from 71 and another from 72, that should agree to within 
the random errors. 

To find the best fit constant bias for this sample we 
perform both the x 2 an d the zero-lag covariance methods 
discussed in sections 14.4.11 and 14.4.21 The curves in the left 
panel of figure [3] show the reduced \ 2 as a function of bias 
for our four smoothing cases. It is worth noting that the 
same calculations for no smoothing does not give a good fit 
(or a minimum) and stays substantially out of the range of 
the plot for any value of the bias. We also show, in the back- 







Bias per slice: 




Smoothing 


z e t f = 0.36 


z eff = 0.70 


z eff =0-62 




[ < z a < 1 ] 


[ K z s < 2 ] 


[ < z s < 2 ] 


Gauss 


1.20 ±0.19 


1.85 ±0.22 


1.70 ±0.17 


TSIS 


0.99 ±0.14 


1.62 ±0.18 


1.47 ±0.13 


NN 


1.59 ±0.22 


2.82 ±0.34 


2.51 ±0.25 


MEF 


0.73 ±0.12 


1.34 ±0.19 


1.19 ±0.13 



Table 1 . Best fit values and errors for each of the four smoothing 
that we have investigated using the covariance method outlined 
in section [4.4.21 Results are shown for the same redshift slices as 
in figure [3] These have been constructed by dividing the source 
(i.e. lensed) galaxies using the photometric redshift (z s ). For each 
of these redshift slices we also show the effective redshift (z e jj). 

ground, shaded regions that come from the measurements 
using the cross-correlations method. We see that the two 
methods for measuring the bias give good agreement. The 
best fit bias does depend on the smoothing scheme that we 
adopt and that in each case we typically measure the bias 
to a precision of 10% to 15%, once a smoothing scheme has 
been chosen. Specific best-fit results and the one sigma er- 
rors for the zero-lag covariance method are shown in table 
[1] We also restate that for the lensed galaxies in the range 
< 2 < 1, the effective redshift, i.e. the mean of the lensing 
efficiency function shown in figurefl] is z e ff = 0.36. In figure 
[4] we show the convergence field at z = 1 coming from our 
four reconstructed densities. We see that all the maps show 
the same broad features, but they vary in the details. Thanks 
to the peaked nature of the TSIS smoothing kernel, we see 
that this reconstruction has more small-scale features than 
the Gaussian and MEF maps. Though the nearest neighbor 
convergence field has a variance that is similar to the other 
maps, we can see from the bottom right panel that the pos- 
itive extremes are larger. This can be attributed to the fact 
that the nearest neighbor scheme has a dynamic smooth- 
ing scale, and we are, therefore, able to resolve the highest 
density peaks, which contain many galaxies. 

5.1.2 Constant Bias with 2 < 2 

Although we restrict our density reconstructions to redshifts 
of less than one, we are still able to use the lensing data at 
higher redshifts. This is because the additional lensing signal 
from the mass in the redshift range 1 < z < 2 will effectively 
behave as an additional noise term. Furthermore, this addi- 
tional noise term will be subdominant to the noise coming 
from the random orientation of galaxies and will not have a 
significant impact on our analysis. The only concern would 
be if there was a strong correlation of the mass distribution 
in these two regimes. However, given the distances involved 
(~ Gpc), this cross-correlation is likely to be subdominant. 

The middle panel of figure [3] shows results using the 
lensed galaxies in the redshift range 1 < 2 < 2 and the 
right panel shows results for the redshift range < 2 < 2. 
Here again, the galaxies are divided into bins with widths 
Az = 0.05, and in each redshift slice the galaxies are di- 
vided into 5 bins in predicted shear. The results are also 
summarised in table [T] For each of the redshift intervals, 
we draw similar conclusions to those from the low redshift 
sample. First, there is good agreement between our two ap- 
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Figure 3. The curves show the reduced x 2 fits for constant bias cases. As we have repeatedly highlighted, it is not a surprise that 
the best fitting bias depends on the specific scheme used to smooth the galaxy density field (bx)- The left panel is for galaxies in the 
range < z < 1. The middle panel is for galaxies in the range 1 < z < 2, and on the right we see the results for all galaxies out to a 
redshift of z = 2. The curves show the reduced \ 2 results for each of our four smoothing cases: (i) Gaussian smoothing; (ii) convolution 
with a Truncated Singular Isothermal Sphere (TSIS); (iii) averaging of the fifth Nearest Neighbours (NN); and (iv) a Multiscale Entropy 
Filter (MEF). To calculate the reduced x 2 , the lensed galaxies have first been divided into bins of width Az = 0.05. Each bin has been 
further divided into 5 bins in predicted shear. The vertical shaded regions show the one sigma measurements coming from the zero-lag 
correlation method outlined in section Pi. 1.21 




proaches for measuring the bias and the best-fit bias de- 
pends on the smoothing scheme. We also see that, by using 
the full sample, we are able to measure the bias to a pre- 
cision of roughly 10%, once a particular scheme is chosen 
for going between the galaxy field and continuous field. The 
precision for a given scheme seems to be at this level regard- 
less of which one is considered. We also see that for all cases 
the measured bias increases with the effective redshift (z e ff) 
being probed. Since this trend is independent of our smooth- 
ing scheme, we can conclude that the bias of our underlying 
sample must increase with redshift. 



5.2 Redshift Evolving Bias 

As we outlined in section 21 it is computationally more sim- 
ple to work with the inverse of the bias, n{z), that is ex- 
panded in a series. Given the expansion shown in equation 
1151 the bias is then given by 



b(z) = 



1 



Mo + Mi2/ 



(25) 



where y — z/1 + z. When expressed in this form, we see 
that we must take care in how we explore possible values of 
Ho and since there is a danger that the bias could be- 
come singular, for instance at /no(l + z) = — H\z or negative. 
To guard against these cases, it is convenient to recast our 
variable such that 



b(z) 



bo 



+ i + fy' 



(26) 



Figure 4. The four panels show the convergence (ft) of the COS- 
MOS field at a redshift of 1 for the four smoothing cases that 
we study (see Table [TJ. Each panel shows the central 1.1 square 
degrees of the COSMOS field. The colour range is the same for 
all panels and corresponds to -0.04 (blue) to O.l(red). 



where / is given by / = /ii/^o- Sensible bounds can now 
be easily placed on / to ensure that the bias is positive and 
non-singular. 
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Figure 5. The one sigma measurements for the bias normaliza- 
tion bo = l/fio an< i the term / that controls the redshift evolution. 
This factor is defined as / = fii/fio and can also be expressed in 
terms of the constant bias / = boW- A constant non-evolving bias 
would have / = 0, which is shown with a dotted line. The results 
shown here are for the cases with fifth nearest neighbour (solid) 
and Gaussian (dashed) smoothing. We show the constraints to 
the low redshifts (red), high redshifts (blue) and the full sample 
out to a redshift of 2 (green). 



Figure [5] shows the one sigma constraints on the two 
parameters / and bo for the nearest neighbor (solid) and 
Gaussian (dashed) smoothing schemes. For each case, we 
show the measurements for the three redshift configurations 
shown in table Q] Once again, we see that the exact normal- 
ization of the bias depends on the smoothing scheme, but 
both cases show evidence for an evolving bias. We also see 
that when only using galaxies in the range < z < 1 we are 
able to measure an overall amplitude, but the constraints on 
the evolution of the bias are very weak. When we include 
the higher redshift sample, we see that the evidence for an 
evolving bias becomes very strong and that the constant bias 
option (dotted line at / = 0) is clearly ruled out. 

In the top panel of figure we show the one sigma 
bounds of the bias as a function of redshift (b(z)) for the 
NN and Gaussian smoothing. In each case we also show the 
best-fit biases. To gain an insight into the significance of 
the bias evolution, we explore the expected bias from a sim- 
ple bias model. For this, we use the bias model presented by 
iMo fc White! l|l996l ). The bottom panel shows the mass limit 
in this model that would correspond to the best-fit biases 
shown in the upper panel. We must be careful in interpret- 
ing these curves, since the bias that we measure is averaged, 
using different smoothing prescriptions, in ways that are not 
accounted for in the analytic model. However, we can con- 
clude that, to be consistent with expectations, there needs 
to be a strong evolution in the tracer with redshift. This is 
not unexpected for a flux limited sample. 



Figure 6. The top panel shows the measured bias as a function 
of redshift coming from the one sigma errors in figure [5] Shown 
are the fifth nearest neighbour (cyan) and the Gaussian smooth- 
ing (orange) results. The bottom panel shows the mass threshold 
that would be necessary in a Mo & White (1996) bias model to 
reproduce the best-fit biases shown in the upper panel. 



6 CONCLUSIONS 

We have created a 3D density reconstruction of the COS- 
MOS field out to a redshift of one using a combination of 
the galaxy clustering and weak lensing measurements. This 
empirical approach means that we do not need to make a 
priori assumptions about the relationship between galaxies 
and dark matter. Instead, we measure this directly while 
creating the optimal density field from galaxy tracers. The 
advantage of relying on a local tracer, such as the density 
of galaxies, rather than an integrated tracer, such as weak 
lensing, is that we are able to produce a density field that 
has a high resolution in redshift. 

We find that if we use only the lensed galaxies that are 
embedded in our density field, i.e. galaxies for which we are 
able to make a full prediction of their lensing signal, then 
we are able to measure a constant non-evolving bias to a 
precision of 10% to 15%, once a smoothing scheme has been 
chosen for going from galaxy points to density field. How- 
ever, we also argue that it is possible to use lensed galaxies 
outside of our density field since the extra lensing signal due 
to the mass at redshifts greater than one will have a very 
weak correlation to the mass in our density. For this rea- 
son, this extra lensing signal will act as an additional source 
of noise that will be sub-dominant to the noise from other 
sources. In this way, we are able to measure the bias to 10% 
by using lensed galaxies out to a redshift of 2. This is en- 
couraging as we look forward to future surveys with areas 
much greater than that covered by COSMOS. A good rule of 
thumb is that parameter errors typically scale as the inverse 
square root of the survey area. In this case, the next gener- 
ation of surveys that will have areas over several hundred to 
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thousands of square degrees can expect to reach sub-percent 
level precision using this method. 

We have looked for evidence of an evolving bias. We 
find that for our density field, which is constructed using 
a flux limited sample, does have a strongly evolving bias. 
This trend remains irrespective of the smoothing scheme 
that we use to go from galaxy points to a smooth density 
distribution. 

The work presented here is our first implementation of 
a method that draws a sharp distinction between local trac- 
ers and integrated measures. In our case, the local tracer is 
the number density of galaxies, which we use to construct 
the matter density field, and the integrated measure is the 
weak lensing signal, which we use to calibrate the free pa- 
rameters in the reconstruction. It is then, in principle, easy 
to extend this method to include further local tracers, such 
as x-ray flux, that can be used to mold the density field and 
extra integrated quantities, such as strong lensing and SZ 
measurements, that can be used to test the accuracy of the 
reconstruction. 
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